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Optical parametric process occurring in a nonlinear planar waveguide can serve as a source of 
light with nonclassical properties. Properties of the generated fields are substantially modified by 
scattering of the nonlinearly interacting fields in a photonic band-gap structure inside the waveg- 
uide. A quantum model of linear operator amplitude corrections to amplitude mean- values provides 
conditions for an efficient squeezed-light generation as well as generation of light with sub-Poissonian 
photon-number statistics. Destructive influence of phase mismatch of the nonlinear interaction can 
fully be compensated using a suitable photonic-band gap structure inside the waveguide. Also an 
increase of signal-to-noise ratio of an incident optical field can be reached in the waveguide. 

PACS numbers: 42.50.Dv Nonclassical states of the electromagnetic field, 42.65.Yj Optical parametric oscil- 
lators and amplifiers 



I. INTRODUCTION 



Properties of linear photonic band-gap structures have 
been an object of intensive investigations in the last sev- 
eral years ^ 0. The most typical characteristics of 
these structures are spatial localization of optical modes 
in confined regions of a given structure and high densi- 
ties of these optical modes. Considering nonlinear ma- 
terials, energy of an optical field is mainly contained in 
these localized modes and so a very strong and thus effi- 
cient nonlinear interaction can occur. For example, sec- 
ond harmonic and sub-harmonic generation in photonic 
band-gap structures has been an object of investigation 
in 0) 13 ■ Phase-matching can be tailored in photonic 
band-gap structures so that fulfilment of phase-matching 
conditions for a given nonlinear process is reached and 
an efficient nonlinear process is guaranteed this way. In 
some cases the overlap of nonlinearly interacting optical 
fields and their mutual spatial phase relations determine 
the strength of nonlinear process and properties of light 
obtained in a nonlinear photonic band-gap structure. 

These properties may also be suitable for the genera- 
tion of light with nonclassical properties (squeezed light, 
light with sub-Poissonian photon- number statistics), as 
has been suggested in Up to now, attention has 

been devoted to the generation of nonclassical light in 
nonlinear photonic band-gap waveguides. It has been 
shown that the process of second-harmonic generation in 
a planar nonlinear waveguide with a corrugation on the 
top can be used to control squeezing of the fundamental 
field 1^ ; the corrugation reproduces a photonic band-gap 
structure. In periodicity of the grating was selected to 



give rise to a longitudinal confinement of the pump field, 
phase matching of the nonlinear process was achieved 
introducing a spatial modulation of nonlinear suscepti- 
bility. Conditions for an efficient squeezed-light gener- 
ation as well as generation of light with sub-Poissonian 
photon-number statistics have been analyzed in for a 
nonlinear waveguide with optical parametric process; the 
photonic band-gap structure was set to assure longitudi- 
nal confinement of the down-converted fields. 

In this contribution we extend the analysis given in 
to account also for the pump-field longitudinal confine- 
ment. This confinement considerably changes amplitude 
and phase relations along the waveguide thus providing 
new possibilities for nonclassical-light generation. Opti- 
cal fields participating in the nonlinear interaction are 
described using the generalized superposition of signal 
and noise. 

A quantum derivation of the equations governing the 
evolution of the interacting optical fields is given in Sec. 
2. In Sec. 3, conditions for squeezed-light generation are 
analyzed. Sec. 4 is devoted to photon-number statistics 
of the generated fields. Possibility to improve signal-to- 
noise ratio of an incident optical field is discussed in Sec. 
5. Sec. 6 contains conclusions. 



II. QUANTUM DESCRIPTION OF THE 
NONLINEARLY INTERACTING FIELDS 

A quantum description of nonlinearly interacting opti- 
cal modes requires the construction of an appropriate mo- 
mentum operator G{z), which then determines Heisen- 
berg equations of motion: 
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FIG. 1: Scheme of the considered nonhnear planar photonic 
band-gap waveguide; App, Asf, and Aip [Aps, Asb, and 
Ais] denote amphtudes of forward-propagating [backward- 
propagating] pump, signal, and idler fields; Kp, Ks, and Ki 
axe linear coupling constants between pump, signal, and idler 
fields, and Kp (Kb) stands for a nonlinear coupling constant 
among forward- (backward-) propagating fields. 



X stands for an arbitrary operator, h is the reduced 
Planck constant, and [ , ] means a commutator. 

If a nonlinear interaction involves counter-propagating 
fields, we cannot straightforwardly assign any momentum 
operator G to the system of interacting optical fields. 
However, we can proceed as follows 8]. We assume the 
nonlinear interaction among all involved fields as if they 
co-propagate and write the momentum operator G{z) in 
the form: 

G{z) = ^ h{'ka)zalaa + ^ h{'ka)zalaa 



a — SB ■'i'E 



\hKs exp{i6iz)al^asB + exp(i(5;z)a|^ai^ 



(2) 



Symbol da {a\) denotes an annihilation (creation) oper- 
ator of mode a; (k(j)z is the corresponding wave-vector 
of mode a along z— axis. We consider six nonlinear ly 
interacting optical modes in the investigated waveguide 
with optical parametric process; forward-propagating 
signal mode (denoted as sp), backward-propagating 
signal mode (ss), forward-propagating idler mode 
{ip), backward-propagating idler mode («_b), forward- 
propagating pump mode (pf), and finally backward- 
propagating pump mode {pb)- Constants Kg, Ki, 
and Kp describe a linear exchange of energy between 
forward- and backward-propagating signal, idler, and 
pump modes. This exchange of energy originates in 
scattering of fields in a photonic band-gap structure. 
Frequency 5i determines periodicity of the corrugation 
on the top of the waveguide. Constants Kp and Kb 
stand for nonlinear coupling coefficients for forward- and 
backward-propagating fields. Values of these parameters 
are determined using simple expressions after a mode 
structure of the waveguide is found (see, e.g., in [gllToj). 
Scheme of the waveguide with the considered interactions 
among the involved optical fields is shown in Fig. ^ 

We now substitute creation operators (d^) of the 
backward-propagating fields by newly introduced auxil- 



iary annihilation operators (6) and vice versa, i.e. 
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(3) 



Heisenberg equations in Eq. Q then have the form: 
dd„^ 



dz 

ddip 
dz 



dz 



db] 



dz 

ddpp. 
dz 



d4^ 
dz 



i(ksj,)zdsp + iKsexp{iSiz)bl^ 
+ 2Kpdppd\^, 

i(ki^)zdjp +iKieyi-p{i5iz)b\^ 
+ 2Kpdppd\^, 

-i{]&SB)zb\B ~ exp(-i(5/z)asj, 
-2KBbl^hB, 

-i(ki^)z6j^ - iK* exp{-iSiz)dip 
-2KBblJsB, 

i{'kpp)zdpp + iKpex.p{2iSiz)blg 
2K pdsp&ip , 

-i(kpB)zbl^ ~ iK*e-Kp{~2i6iz)dp 



2K*Bbl,b 



(4) 



To reach the final equations, we have to make the fol- 
lowing steps: 1. Return to the original operators d^, 
d in Eqs. I0J using the substitution in Eq. 2. 
Transform Eqs. Q into the interaction picture {da{z) = 
Aa{z) exp[i{'ka)zz])- 3. Write operators Aa{z) in the 
interaction picture as Aa{z) — Aa{z) + dAa{z), where 
Aa{z) is a classical amplitude mean- value and dA{z) is a 
small operator correction to this amplitude mean-value. 
This procedure results in a system of nonlinear differen- 
tial equations for amplitude mean- values Aa and a system 
of linear operator differential equations for small operator 
amplitude corrections SAa- 

The system of nonlinear differential equations for clas- 
sical amplitude mean- values Aa is written as follows: 



dA, 



dz 

dA^ 

dz 

dAsB 
dz 

dAi„ 



dz 



iKs eyip{-i5sz)AsB 

+ 2Kp eKp{i5pz)AppA*^, 

iKi exp{-i5iz)AiB 

+ 2Kp eKp{i5pz)AppA*^, 

-iK* exp{i6sz)Asp 

- 2KBexp{~idBz)ApBA*^, 

-iK* eyip{i6iz)Aij, 

- 2KBexp{-iSBz)ApBAl^, 
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dAr. 



dz 

dA. 



PB 



dz 



iKp exp(— i(5pz) 

- 2Kpe-Kp{~iSFz)AspAip, 
—iK* exp(i(5pz)ylp^ 

+ 2K*g exp(i(5sz)yls^Ajg , 



(5) 



and 



l(k.J.|-|(k.J 



b^F,B. (6) 



Any solution of the system in Eqs. |SJ| obeys the fol- 
lowing conservation law of energy (in quantum interpre- 
tation, "the overall number of virtual photons in the in- 
teraction" is conserved): 



-1^ 



\A, 



2\A 



PB I 



= 0. 



(7) 



If the nonlinear terms in Eqs. (|SJ) are omitted, the 
solution of Eqs. (jSJ can be written as: 



4(0) 



exp —I 



.SaZ 



Ba cos{Aaz) + Ba sin(Aa2:^ 



exp I 



.6aZ 



^1 ( ^ttS^ cos{Aaz) + sin(Aaz; 
sin(Aa2;) - cos{Aaz] 



2Ka 



a — s,i,p 



and 



A. 



(52 

-f-\Ka\'; 



a — s, I, p. 



(8) 



(9) 



In Eqs. constants Bg, Bg, Bi, B^, Bp, and Bp are 
set according to boundary conditions at both sides of the 
structure. The solution of the nonlinear set of equations 
written in Eqs. © is reached numerically using an it- 
eration from the solution written in Eqs. ijSJi. A finite 
difference method called BVP has been found to be 
suitable for this task. 

The evolution of small operator amplitude corrections 
SAa is governed by the following equations: 



dSA 



dz 
d6A, 



dz 
dSA, 



= ICsSAsB 

= JCiSAig 



dz 



IClSAs 



- ICf 

- ICb 



ApJAl+AlSAp, 



AppSAl^ + A^^SApj, 



Ar)„6A,: 



A*JA„ 



dSA, 
dz 
dSA^ 



IC*5A,, 



^ - K,p6Ap, 



d6A. 



JCb 



ApB^Al^ + A*^SAp^ 
AspSAi^ ~\- AipSAsf, 
As a SAi „ -|- Ai „ 6Af: „ 



(10) 



Functions ICs, JCi, ICp, ICf, and ICb introduced in Eqs. 
((Tn)l are defined as: 



ICa = iKaexp{-iSaz), 
K-F = 2Xi? exp(i(5i?z), 
ICb = 2Kb CTq){—iSBz). 



a = s,i,p, 



(11) 



The solution of the system of linear equations in Eqs. 
(|10|l for operator amplitude corrections SAa can be found 
numerically and written in the following matrix form: 



SAf^ouI 

SAbm 



Uff 
Ubf 



Ufb 
Ubb 



f SAfm 
\ 5AB,oMt 



(12) 



where 



SAfm = 



5Ai 



/'5i.,(0)\ 




^SAsAL)\ 


5AIA^) 




SAUL) 




, SAF,out = 


SAAL) 

SAlAL) 


Mp,(0) 




SApAL) 


\5AIAQ)) 




\SAIAL)J 


fSAsjL)\ 




fSAsAO)\ 


SAIJL) 




sAUm 


SA.AL) 
SA^^AL) 


, SAb, out = 


5 A, AO) 
SAliO) 


SApAL) 




SApAO) 


\SAI{L)J 




\SAl{0)J 



(13) 



Matrices Uff, Ufb, Ubf, and Ubb characterize the so- 
lution of Eqs. l(Tn)l. 

Input-output relations among linear operator ampli- 
tude corrections SAa can be found solving Eqs. H12(l with 
respect to vectors SAF.out and SAB.out- 



SAF,out 



/ Uff — UfbUqqUbf UfbU^q 



U 



UurUbF 



BB 
5AF;in 
5AB:in 

5Af., 
5Ab 



U 



BB 



(14) 
(15) 



Matrix U defined in Eq. (|15|l describes input-output re- 
lations among the linear operator amplitude corrections 
SAa- The output linear operator amplitude corrections 
contained in vectors 5AF,o\it and 6AB,out obey boson 
commutation relations provided that the input linear op- 
erator amplitude corrections occurring in vectors 5AF,in 
and SAbaii obey boson commutation relations. It has 
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been shown in 1^ that this nontrivial property is fulfilled 
by any system described by a quadratic hamiltonian. 

The method of derivation of the operator equations in 
Eqs. (|10|l through the set of operator equations written 
in Eqs. I@J reveals that the following "commutation re- 
lations" among the small operator amplitude corrections 
SAaiL) are fulfilled: 



[SA,{L),SAk{L)] = 


0, 


[SML),SAl{L)] = 


^ifc, 


[dML),6Al{L)] = 


0, 


[5Ai{L),6A-,{L)] = 


0, 


[SA,iL),SAliL)] = 




[SAj{L),SA-,{L)] = 


0, 


i,k = SF,iF,PF, 


i, k 



(16) 



These relations have been found to be useful in control- 
ling precision of the numerical solution. 

We describe the interacting fields in the framework of 
the generalized superposition of signal and noise ^3 (co- 
herent states, squeezed states as well as noise can be con- 
sidered). Any state of a two- mode field is determined by 
values of parameters Bj, Cj, Djk, and Djk 



D 
D 



{AA]AA,), 

((Ai,r), 
(Ai.Aifc), 



(17) 



AAj = Aj — (Aj). Symbol ( ) stands for a quantum 
statistical mean value. Coefficients Bj, Cj, Djk, and Djk 
can then be determined using matrix U introduced in Eq. 
()15|l and incident values of i3j,in,^ and Cj^n.yt related to 
anti-normal ordering of field operators (for details, see 



Cj,in,A = ^exp(z?9j)sinh(2rj). 



(18) 



Symbol rj denotes a squeeze parameter of the incident j- 
th mode, 'dj means a squeeze phase, and rich.j stands for 
a mean number of incident chaotic photons. Coefficients 
Djk, in, A and Djk, in, A for an incident field are set to zero 
because the incident fields are assumed to be statistically 
independent. 

The expressions for coefficients Bj, Cj, Djk, and Djk 
can be written in terms of matrix elements of Z// as follows 
i: 

6 

Dj ^ [(^27-1.2fc-1^2j~l,2fcC'fc -I- C.C.) 

k=l 

+ |^2j-l,2fc-lP {Dk,in,A ~ 1) + \i^2j-l,2k\'^ Bk^in,A 



u 

Cj = ^2j-l,2k-l^k,in,A + ^2i-l,2k^k,in,A 

fc=l 

+ 1^2 j-l,2k-ll^2j -1.2k {2BkAn,A — 1) , 
6 

^J*; ~ [^^2j-l,2i-1^2fc-l,2i-lC';,in,.4 
1=1 

+ W2j-l,2;W2fe-l,2/C';*in,X 
+ l^2j-l,2l-ll^2k-l,2lBlAn,A 

+ U2]-l,2lli2k-l,2l-l{BlSn,A- , 3 ^ k, 
6 

\^-l^2j-l,2l^2k-l,2l-lCl,in,A 



Djk 



1=1 



' ^2j-l,2i-1^2fe-l,2(C';*in,^ 

' ^2j-l,2i-1^2fe-l,2i-l (Si,in,^ " 1) 



^2i-l,2l^2k-l,2lBlin,A 



, 3 J^k; 



(19) 



C.C. stands for complex conjugated terms. 



III. SQUEEZED-LIGHT GENERATION 

The level of noise present in quadrature components qj 
[(jj = Aj + Aj , Aj stands for an electric-field-amplitude 

operator of mode j] and pj [pj — —i{Aj — A'j)] appro- 
priate for mode j can be lower than the level of noise 
characterizing the vacuum field. Then we speak about 
squeezed light. In general, the maximum amount of 
available squeezing is reached under some chosen value 
of a local-oscillator phase in the homodyne-measurement 
scheme and the corresponding amount of squeezing is 
given in theory by principal squeeze variance Xj |14| . 

It is useful to combine some optical fields on a beam- 
splitter and to study properties of the output fields. Such 
fields can have a nonclassical character under certain con- 
ditions. For example, signal and idler fields generated in 
optical parametric process have this property, because 
one signal photon and one idler photon are created to- 
gether in one elementary quantum event of the nonlinear 
process. We use the notation compound mode for this 
case and define the appropriate operators for quadrature 
components combining j-th and A:-th modes; qjk = (jj+qk 
and pjk = Pj +Pk- 

Using the parameters characterizing a state and de- 
termined in Eqs. (|19|l . we obtain for single- mode prin- 
cipal squeeze variance Xj and compound-mode_^rincipal 
squeeze variance Xij the following expressions 



A, 



1 + 2[Bj 
2\1 + B 



\C 



Cj\h (20) 
J , Sfc - 2Re{D,k) 
Ck + 2Djk\]. (21) 



Values of principal squeeze variance Xj less than one in- 
dicate squeezing in a single-mode case. Squeezed light is 
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generated in a compound-mode (two-mode) case if values 
of principal squeeze variance Aj^ are less than two. 

Similarly as in 0, assuming KgL, KiL, KpL, KpAL 
{A being a typical field amplitude), and KbAL being 
small, analytical expressions for principal squeeze vari- 
ances can be found solving equations in Eqs. ((Tn|l iter- 
atively. The obtained expressions for principal squeeze 
variances for single-mode and compound-mode cases are 
the same as those given in Eqs. (23) and (25) of 0, 
where the system with Kp = is analyzed. Thus, no in- 
formation about the influence of linear pump scattering 
on squeezed-light generation can be obtained. 

In the following discussion a strong incident forward- 
propagating pump field and also nonzero incident 
forward-propagating signal and idler fields are consid- 
ered. Squeezed light cannot be generated in single-mode 
cases. However, in general, compound modes {sp,ip), 
(sb,*b), and {sFjis) generate squeezed light under cer- 
tain conditions. We note, that properties of the optical 
fields that follow from the symmetry between the sig- 
nal and idler fields are not mentioned explicitly. For 
example, following the symmetry, also compound mode 
(«i?,ss) can be squeezed. It is interesting to note, that 
the analyzed waveguide behaves qualitatively similarly as 
a waveguide with two separated parts with optical para- 
metric processes and with fields in different parts inter- 
acting linearly through evanescent waves. This waveg- 
uide with co-propagating fields has been analyzed in 

Cllifij. 

We first consider linear scattering only between the 
pump modes {Kp 7^ 0). If there is phase matching of all 
interactions, the greater the value of Kp the greater the 
values of principal squeeze variance A of mode (sFjip)- 
This means that phases of the nonlinearly interacting op- 
tical fields along the structure are modified by nonzero 
values of Kp in such a way that the amount of gen- 
erated squeezing decreases. On the other hand, linear 
coupling given by nonzero values of Kp transfers energy 
into mode pb and so squeezing can be observed also in 
mode (sB,iB). However, greater values of Kp suppress 
the amount of squeezing in this mode, as is shown in Fig. 
11 

In general, the greater the pump-field amplitude 
the lower the value of principal squeeze variance A. Non- 
linear coupling constants Kp and Kb behave in the same 
way. Also the dependence of principal squeeze variances 
A on length L of the waveguide is similar, because "the 
overall amount of nonlinear interaction" is proportional 
to L. 

If the waveguide is characterized by a greater value of 
Kp then the value of A depends strongly on linear phase 
mismatch Sp of the forward- and backward-propagating 
pump fields. A dramatic decrease of principal squeeze 
variance A in mode (sfj^f) for 6p > '2\Kp\ is shown in 
Fig. 121 As the approximate solution in Eqs. ((SJ for am- 
plitude mean-values suggests, this region of parameters 
with low values of A can be characterized by an oscillating 
behavior of the amplitudes (as functions of the position 




FIG. 2: Principal squeeze variance A of mode {33,13) as a 
function of forward-propagating pump amplitude App and 
linear coupling constant Kp between the pump modes; Kp = 
Kb = 5 X 10-2, Ks ^ Ki = 0, 5» = (5, = 0, Sp = 0, 
Sp = Sb = 0, Asjr = = 0.1, Apg = Asi, = Aig = 0; 
incident coherent states are assumed. The following units 
are used for the considered physical quantities: [Kp] = 
[Kb] = 10-«mm-imV-\ [Ks] = [K] = [Kp\ = mm'S [L] = 



mm, [Ss] = 



= [5v] = [5f] = [Sb] = mm-\[Aa] 



^ 1.0 



0.5 



20 40 60 BO 100 

5„ 



FIG. 3: Principal squeeze variance A of mode {sp,ip) as a 
function of linear phase mismatch Sp between the pump fields; 
L — 2, Kp = 5, App = 10, and values of the other parameters 
are the same as in Fig. 2. 



in the waveguide) . A similar decrease of A occurs also in 
mode {sbtIb) as 5p exceeds 2|ifp|, but for greater values 
of 6p an increase of values of A follows because transfer 
of energy to the backward-propagating pump field de- 
creases. 

If the nonlinear interaction is not phase-matched 
{5ni ^5p = 5b^Q), then the linear scattering of the 
pump field (described by Kp and 5p) can compensate for 
nonzero values of nonlinear phase-mismatch 5ni and en- 
able better values of principal squeeze variances A. If, 
for example, we consider 5ni — 5 mm"^ and values of 
the other parameters written in caption to Fig. 21 prin- 
cipal squeeze variance A of mode {sp,ip) is greater than 
0.8 assuming the waveguide without a photonic band-gap 
structure. Having a photonic band-gap structure with 
suitable values of parameters inside the waveguide, prin- 
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FIG. 4; Topological graph of principal squeeze variance A in 
mode (sFjip) as it depends on linear coupling constant Kp 
and linear phase mismatch Sp between the pump fields; L = 2, 
App — 10, Sf ~ Sb = 5, and values of the other parameters 
are the same as in Fig. 2. 



FIG. 5: Principal squeeze variance A of mode (spjip) as a 
function of linear phase mismatches Ss and Si between the 
signal and idler fields, respectively; L — 2, Kp = 5, Ks — 
Ki — 5, App — 10, and values of the other parameters are the 
same as in Fig. 2. 



cipal squeeze variance A can even reach the value 0.04 
appropriate also for a phase-matched nonlinear interac- 
tion in the waveguide without a photonic band-gap struc- 
ture. Suitable values of Kp and Sp lie typically around 
the curve 6p = 2\Kp\, as can be seen in Fig. ^ Also 
greater values of Kp are needed to have lower values of 
A. The possibility to compensate for a nonlinear phase- 
mismatch using linear coupling between modes has been 
discussed in |l7l | from the point of view of efficiency of 
energy conversion. 

Nonlinear phase mismatch 6ni can be also compen- 
sated using nonzero incident amplitudes of the backward- 
propagating pump mode. We have found the opti- 
mum value of the phase difference between the forward- 
propagating pump amplitude and backward-propagating 
pump amplitude to be Stt/2 [tt/2] for mode (sfj^f) 
[(ss,«b)] with respect to squeezed-light generation. 

Now we assume linear scattering of pump, signal, and 
idler fields {Kp ^ 0, K^ 0, K, =^ 0). The analysis of 
squeezed-light generation is difficult under these condi- 
tions because every linear coupling constant introduces 
some phase changes along the waveguide that modify the 
nonlinear interaction. In general, squeezing can be ob- 
served in modes {sF^ip), (sb,?^), and (sf,«b) under 
suitably chosen values of parameters. The occurrence of 
squeezing in mode (sf, is) originates in linear scattering 
of the down-converted fields. 

Nonzero values of Kg and Ki support squeezed-light 
generation in mode {33,13), because they transfer energy 
from modes sp and ip into modes S3 and 13 and so the 
nonlinear interaction among the backward-propagating 
fields can have also a strong stimulated part. A nonzero 
value of Ki is also indispensable for observing squeezing 
in mode {sp, 13), because "an already generated squeezed 
light in the nonlinear interaction among the forward- 
propagating fields has to be transferred into mode 13" . 

Being in non-oscillating regime of behavior of ampli- 
tude mean-values {5p < 2\Kp\, Sg < 2\Ks\, Si < 21X^1) 



lower values of A can be reached in mode {33,13) in com- 
parison with those reached in mode {3F,ip) probably be- 
cause modes 33 and is begin the nonlinear interaction 
in the vacuum states and so phases of their amplitudes 
can be suitably set in the interaction. In general, the 
regime with the oscillating behavior of amplitude mean- 
values {Sp > 2\Kp\, Ss > 2\Ks\, S, > 2\Ki\) is better 
for squeezed-light generation and low values of princi- 
pal squeeze variances A can be reached in this regime. 
Squeezing can be observed under a wide range of values 
of parameters in non-oscillating regime in modes {sp, ip), 
{33,13), and (sFjis). We note, that greater values of, 
e.g.. Si effectively decrease the value of Ki which acts 
against squeezed-light generation in mode (sfj^b)- A 
difference in the ability to generate squeezed-light in 
oscillating and non-oscillating regimes is clearly visible 
in Fig. |S1 where principal squeeze variance A of mode 
(sf,*s) is plotted as a function of linear phase mis- 
matches Sg and Si. 

As can be shown directly by a suitable substitution 
into Eqs. l(SJl and H10|) . their solutions depend only on 
the overall phase ■!/; = — arg(_K'p) -I- aig{Ks) + a,ig{Ki) of 
the linear coupling constants. In our investigations, prin- 
cipal squeeze variances A in modes {3F,ip) and {33,13) 
had minimum values for i/) = 7r/2 and A in mode {3F,i3) 
reached its lowest value for t/j = 37r/2. Also the depen- 
dence of principal squeeze variances A on phases of the 
incident forward-propagating signal and idler fields have 
not been observed. 

At the end of the discussion of squeezing we men- 
tion the case in which, e.g., the incident forward- 
propagating signal field is stronger than the incident 
forward-propagating pump field. We then have squeez- 
ing also in compound modes that combine one pump field 
with one down-converted field under suitably chosen val- 
ues of parameters. In general, values of principal squeeze 
variances A are lower in these compound modes that con- 
tain a down-converted field with great values of ampli- 
tudes inside the structure. 
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IV. PHOTON-NUMBER STATISTICS 

When a field is detected by a classical detector statisti- 
cal properties of the incident field are suitably described 
by normally-ordered moments of integrated intensity W 



2,3, 



(22) 



where A denotes an electric-field-amplitude operator of 
the incident optical field. 

Type of a statistical distribution of photoelectrons 
emitted inside the detector can be determined using Fano 
factor Fn- Fano factor Fn can be expressed in terms of 
normally-ordered moments of integrated intensity W as 
follows: 



F„, 



((An)^) 
in) 



1 



(23) 



Symbol n denotes the number of photoelectrons. An — 
n — (n), and AW = W — {W)j^. Intensity operator Wij 
of a compound mode is then determined along the 
relation W^j = + Wj, where Wi (Wj) stands for the 
intensity operator of mode i (j). Classical fields obey the 
inequality Fn > 1. On the other hand values of Fano 
factor Fn smaller than one can be reached considering 
nonclassical fields (sub-Poissonian light). The condition 
Fn < 1 means that fluctuations in the number of pho- 
toelectrons are suppressed below the classical limit that 
is given by the Poissonian photon-number statistics of a 
laser radiation. 

Assuming that the outgoing fields can be described in 
the framework of the generalized superposition of signal 
and noise, moments of integrated intensity W can be 
written as: 



2 



c.c. 



(24) 



Coefficients Bj, C,-, Djk, and Djk are given in Eqs. (|19|l . 
Symbol S^j in Eqs. (pljl denotes a coherent amplitude 
of the jth outgoing field. We note that photon-number 
distribution as well as moments of integrated intensity 
can be determined in general usin g an expansion into 
Laguerre polynomials (see, e.g., [l3lll8l|). 

Nonclassical character of photon-number statistics oc- 
curs mostly at single-photon level. For this reason, we 
assume a regime in which the waveguide is pumped by 
a strong forward-propagating pump field and classical 
strong amplitude mean- values Agp , Ai^ , Agg , and Ai^ 
are zero. Operator amplitudes A of the signal and idler 
fields are then given just by their linear operator ampli- 
tude corrections 5A. We also assume that the incident 
fields described only by linear operator amplitude correc- 
tions are coherent and denote their amplitudes by ^. 




FIG. 6: Fano factor of mode {sp,ip) as a function of 
forward-propagating pump amplitude App-, input linear op- 
erator amplitude corrections 5 A are assumed to be in coherent 
states with amplitudes ^, amplitudes ^ are expressed in units 
of 10 Vm~^ and so mean values of intensities are directly equal 
to mean photon numbers; L = 2, Kf = Kb = 5 x 10"'^, 
Kp ^ Ks = = 0, 5p = Ss = = 0, Sf ^ Sb = 0, 
App = 10, Asp = Aip = 0, Apg = Asg = Aig = 0; 

CsF — ~10, ^ip ~ 10 ^pp ~ ~ £,iB ~ CpB ~ 0- 



Assuming KgL, KiL, KpL, KpAL {A being a typical 
field amplitude), and KbAL being small, expressions for 
first and second moments of integrated intensity W can 
be found as it was done in ■ This approximation gives 
the same expressions as those published in Eqs. (30) and 
(31) in 01 for the case Kp = 0; i.e. no conclusion about 
the influence of Kp can be deduced. 

We first consider a phase-matched nonlinear waveg- 
uide with no photonic band-gap structure {K^ = Ki = 
Kp = 0). Sub-Poissonian light in mode (sFjip) can 
occur for a sufficiently strong pumping. However, as 
Fig. shows, if the value of pump amplitude App is too 
great, sub-Poissonian character of the generated light is 
lost. Suitable values of pump amplitude App for sub- 
Poissonian-light generation depend on length L of the 
waveguide. The longer the waveguide, the smaller the 
suitable values of pump amplitude App . A photonic 
band-gap structure with Kp ^ (also Ks = Ki = 
is assumed) inside the waveguide in this case leads to 
a redistribution of energy in the pump modes in the 
way that enables again sub-Poissonian-light generation 
in mode (sf,*f)- For a given value of length L and 
a given sufficiently great value of App there is an op- 
timum value of Kp for which the value of Fano factor 
Fn reaches a minimum value that is obtained also in 
the waveguide without a photonic band-gap structure 
{Fn ~ 0.3, see Fig. [T)). If the value of pump ampli- 
tude App is small, the greater the Kp the greater the 
Fano factor F„; i.e. a photonic band-gap structure does 
not support sub-Poissonian-light generation in this case. 
Nonzero values of linear phase mismatch Sp in this oth- 
erwise phase-matched interaction result in greater values 
of Fano factor Fn ■ 

In order to obtain sub-Poissonian light in mode 
(sfi^f), the nonlinear interaction has to be stimulated, 
i.e. the incident small amplitudes and have to 
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FIG. 7: Fano factor F„ of mode {sp,ip) as a function linear 
coupling constant Kp between the pump fields; App = 10 and 
values of the other parameters are the same as in Fig. 6. 
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FIG. 8: Fano factor Fn of mode (sf,*f) as a function of 
amplitude ^: = ^ip = Kp = 1.4, App = 10, and 
values of the other parameters are the same as in Fig. 6. 



be nonzero. Moreover, values of amplitudes ^^'^ Cif 
have to be nearly the same and their phases have to fulfill 
the condition arg(^SF) + arg(^ip) w tt. Under these con- 
ditions Fano factor F„ of mode (sFjip) decreases with 
increasing values of amplitudes S,sp and S,ip . For a certain 
value of these amplitudes a minimum value of Fano factor 
Fn is reached. In our case this occurs for approximately 
100 incident photons in both forward-propagating signal 
and idler fields, as is shown in Fig. |H1 

If the nonlinear interaction is not phase-matched 
(dni = Sp = Sb ^ 0) in the waveguide with no photonic 
band-gap structure {Kg = Ki = Kp = 0) sub-Poissonian 
light in mode {sF,iF) can still be obtained but greater 
values of Fano factor F„ occur. Even values of Fano fac- 
tor Fn can monotonically increase as the value of pump 
amplitude App increases in some cases. Then values of 
parameters of the photonic band-gap structure (Kp ^ 0, 
6p 7^ 0, Kg = Ki = 0) can be set in such a way that 
the original low values of Fano factor F„ are restored. If 
we consider Sni = 5 mm^^ under the conditions specified 
in caption to Fig. the lowest value of Fano factor F„ 
is approximately 0.8. A suitable choice of values of Kp 
and Sp provides the original value of Fano factor Fn be- 
ing roughly 0.3, as is documented in Figl^l According to 
our investigations, the region in {Kp, Sp) space for which 
the influence of nonlinear phase mismatch Sni is compen- 
sated lies around Sp « 2\Kp\ and also greater values of 
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FIG. 9: Topological graph of Fano factor F„ of mode (sFjip) 
in dependence on linear coupling constant Kp and linear phase 
mismatch Sp between the pump fields; only values of F„ lower 
than 1 are plotted; L = 1, Sp = Sb = 5, App = 10, and values 
of the other parameters are the same as in Fig. 6. 




FIG. 10: Fano factor F„ of mode (sf,*f) in dependence on 
linear coupling constants Ka and Ki between the signal and 
idler fields; Kp = 1.4, App = 10, and values of the other 
parameters are the same as in Fig. 6. 

Kp are needed. As is seen from Fig. O the region where 
Sp is greater than 2|i4rp| is not suitable for the generation 
of sub-Poissonian light. 

A suitable choice of the incident phase ippp of the 
forward-propagating pump mode [(ppp — a.rg{App)] can 
be used as a final step to reach the lowest possible value 
of Fano factor F„ that is allowed by a given set of val- 
ues of waveguide parameters. We have observed a strong 
dependence of values of Fano factor F„ on the value of 
phase ippp . 

Low values of Fano factor Fn are usually reached when 
also integrated intensity W oi a, given mode has low val- 
ues. 

We now consider the influence of photonic band-gap 
structure in its full complexity, i.e. also the signal and 
idler fields are linearly scattered {Kg ^ 0, Ki ^ 0). 
Under phase-matched conditions increasing values of Kg 
and Ki destroy sub-Poissonian photon-number statistics 
in mode {sftIf) (see Fig. [TU)) . On the other hand, 
they support sub-Poissonian-light generation in modes 
{sB,iB) and {sF,iB)- Mode {sF,iB) can have values of 
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Fano factor F„ less than one only for small values of 
Kg and also a nonzero value of Ki is required. In gen- 
eral, oscillating regime of behavior of field amplitudes 
{5s > 2\Ks\, Si > 2\Ki\, Sp > ^IKpl) is more convenient 
for sub-Poisonian-light generation. 

Values of Fano factors Fn in general depend only on 
the phase ip [ip — — arg{Kp) + arg{Ks) + aTg{Ki)] that 
combines phases of all linear coupling constants. The 
lowest value of Fano factor Fn in mode (s^, is) is reached 
for tjj = 7r/2; ijj = 37r/2 has been found to be optimum 
for mode (sf, is)- 

We note that a qualitative similarity can be found in 
the behavior of photon-number statistics between the in- 
vestigated waveguide and that one composed of two sep- 
arated nonlinear parts with co- propagating fields [l^ . 



V. INCREASE OF INTENSITY 
SIGNAL-TO-NOISE RATIO 




FIG. 11: Topological graph of second reduced moment Rw 
of integrated intensity in mode sf as a function of linear cou- 
pling constant Kp and linear phase mismatch 5p between the 
pump fields; an incident value of Rw equals 1.75; App = 10, 
Sf ^ Sb ^ 5, = C^F = 10, Uch.sp = 100, nch,ip = 0, and 
values of the other parameters are the same as in Fig. 6. 



The nonlinear waveguide can also be used for improv- 
ing signal-to-noise ratio of an incident field. This effect 
can be explained claiming that the signal part and the 
noisy part of the incident field have different amplifica- 
tion coefficients in the nonlinear process. An average 
amplification coefficient of a field in a nonlinear pro- 
cess depends on a statistical distribution of incident-field 
phases; there exists one phase for which the amplifica- 
tion is maximum. If the central phase of the incident 
field (corresponding to a coherent signal amplitude) has 
the strongest amplification then the noisy part (with a 
blurred phase distribution) is less amplified on average 
and signal-to- noise ratio increases. 

Similarly as for sub-Poissonian-light generation, the 
best conditions for the reduction of noise-to-signal ra- 
tio can be found in a phase-matched waveguide with- 
out a photonic band-gap structure. If nonlinear phase- 
matching cannot be reached, a suitable photonic band- 
gap structure inside the waveguide compensates for non- 
linear phase-mismatch and gives similar conditions for 
the reduction of noise-to-signal ratio as those found in the 
perfectly phase-matched waveguide. For a given value of 
nonlinear phase mismatch Sni (Sni = 6p = Sb), there are 
regions in space {Kp,5p) where the optimum conditions 
for the reduction of noise-to-signal ratio are found. These 
regions lie around the line 5p = '2\Kp\ and also values of 
Kp have to be greater. Second reduced moment of inte- 
grated intensity Rw [Rw = (W"^) I for mode sp 
having incident Rw = 1.75 (100 incident signal photons, 
100 incident noisy photons) is plotted in Fig. Illlassum- 
ing 6ni — 5 mm^^. Regions suitable for the reduction 
of noisc-to-signal ratio are visible in Fig. ^2 the lowest 
achieved value of Rw is 1.36. 

Capability to reduce noise-to-signal ratio increases 
with the increasing forward- propagating pump amplitude 
App (see Fig. I12II . This monotonous behavior clearly 
shows that the nonlinear interaction is responsible for 
this effect. The strength of the effect also depends on 




FIG. 12: Relative second reduced moment Tw of integrated 
intensity in mode sf as a function of forward-propagating 
pump amplitude App and number Uch.sp of incident chaotic 
photons in this mode; Tw = Rw /^w where R}^ [Rw^) 
characterizes the incident (outgoing) field; Kp = 8, (5p = 20, 
5f ~ 5b ~ 5, £,sp — iip — 10, Wch.ip — 0, and values of the 
other parameters are the same as in Fig. 6. 



the amount of incident noise. As demonstrated in Fig. 
1121 for mode SF^ an optimum value of the number of in- 
cident noisy photons exists for which the reduction of 
noise-to-signal ratio is the best. 



VI. CONCLUSIONS 

A planar nonlinear photonic band-gap waveguide with 
optical parametric process has been analyzed from the 
point of view of generation of squeezed light and light 
with sub-Poissonian photon- number statistics. It has 
been shown that squeezed light as well as sub-Poissonian 
light can be generated in compound modes composed 
of one signal and one idler field; both forward- and 
backward-propagating fields can be successfully com- 
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bined. The waveguide with a photonic band-gap struc- 
ture cannot provide better values of principal squeeze 
variance and Fano factor in comparison with the waveg- 
uide with no photonic band-gap structure and having 
the nonlinear interaction phase- matched. However, if the 
nonlinear interaction in a waveguide cannot be phase- 
matched for some reason, inclusion of a photonic band- 
gap structure such as to compensate for phase mismatch 
leads to values of principal squeeze variances and Fano 
factors that are found assuming perfect phase-matching. 
This property makes nonlinear photonic band-gap waveg- 
uides promising as sources of light with nonclassical prop- 
erties. Moreover, if a noisy light is incident on the waveg- 
uide its signal-to-noise ratio can be improved as the light 



propagates in the waveguide. 
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